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ABSTRACT 

We present detailed calculations of the prompt spectrum of gamma-ray bursts 
(GRBs) predicted within the fireball model framework, where emission is due to 
internal shocks in an expanding relativistic wind. Our time dependent numerical 
model describes cyclo-synchrotron emission and absorption, inverse and direct 
Compton scattering, and pair production and annihilation (including the evo- 
lution of high energy electro-magnetic cascades). It allows, in particular, a self- 
consistent calculation of the energy distribution of pairs produced by photon 
annihilation, and hence a calculation of the spectra resulting when the scattering 
optical depth due to pairs, t±, is high. We show that emission peaks at ~ 1 MeV 
for moderate to large r±, reaching r± ~ 10^. In this regime of large compactness 
we find that (i) A large fraction of shock energy can escape as radiation even for 
large r±; (ii) The spectrum depends only weakly on the magnetic field energy 
fraction; (iii) The spectrum is hard, e^dN/de oc with 0.5 < a < 1, between the 
self absorption {sssa = 10°-^±°-5 keV ) and peak (speak = 10°-^±°-^ MeV ) photon 
energy, (iv) and shows a sharp cutoff at ~ 10 MeV; (v) Thermal Comptonization 
leads to emission peaking at Epeak ^ 30 MeV, and can not therefore account for 
observed GRB spectra. For small compactness, spectra extend to > 10 GeV with 
flux detectable by GLAST, and the spectrum at low energy depends on the mag- 
netic field energy fraction. Comparison of the flux at ~ 1 GeV and ~ 100 keV 
may therefore allow to determine the magnetic field strength. For both small 
and large compactness, the spectra depend only weakly on the spectral index of 
the energy distribution of accelerated electrons. 

Subject headings: gamma rays: bursts and theory — radiation mechanisms: 
non-thermal — methods: data numerical and analj^ical 
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1. Introduction 



It is widely accepted that 7-ray bursts are produced by the dissipation of kinetic energy 
in a highly relativistic wind, driven by gravitational collapse of a (few) solar mass object into 
a neutron star or a black hole (see e.g. Piran (2000); Meszaros (2002); Waxman (2003) 
for reviews). The prompt 7-ray emission is believed to be produced by synchrotron and 
inverse- Compton emission of electrons accelerated to relativistic energy by internal shocks 
within the expanding wind (see, however, Lazzati, Ghisellini, Celotti & Rees 2000; Ghisellini, 
Lazzati, Celotti & Rees 2000). Synchrotron emission is favored if the fireball is required to 
be "radiatively efficient," i.e. if a significant fraction of the fireball energy is required to be 
converted to 7-rays. 

Over a wide range of model parameters, a large number of pairs are produced in 
internal coUisions, due to annihilation of high energy photons (e.g. Meszaros & Rees 2000; 
Guetta, Spada & Waxman 2001). In fact, if the internal shocks occur at small enough radii, 
the plasma becomes optically thick and a second photosphere is formed (Meszaros & Rees 
2000; Meszaros, Ramirez- Ruiz, Rees & Zhang 2002) beyond the photosphere associated 
with the electrons initially present in the fireball. As we show here (§ 2.2, see also Guetta, 
Spada & Waxman 2001) requiring the emission to be dominated by ~ 1 MeV photons 
implies, within the fireball model framework, a moderate to large optical depth due to 
scattering by pairs. When the scattering optical depth due to pairs is high, calculation of the 
emergent spectrum becomes complicated. Relativistic pairs cool rapidly to mildly-relativistic 
energy, where their energy distribution is determined by a balance between emission and 
absorption of radiation. The emergent spectrum, which is affected by scattering off the 
pairs population, depends strongly on the pairs energy distribution, and in particular on 
the "effective temperature" which characterizes the low-end of the energy distribution. The 
pairs energy distribution is difficult to calculate analytically. Moreover, analytic calculation 
of the spectrum emerging from the electromagnetic cascades initiated by photon annihilation 
is also difficult. Therefore, in order to derive the emergent spectrum a numerical model is 
needed (see e.g. Ghisellini & Celloti 1999; Zhang & Meszaros 2003). 

Emission from steady plasma, where pair production and annihilation are taken into 
consideration, was numerically studied in the past in the context of active galactic nuclei 
(AGN's). It was found that thermal plasma, optically thin to Thomson scattering, and 
characterized by a comoving compactness 10 < /' < 10^, has a normalized pair temperature 
e = kT/meC^ ^ 10-2 - 10-^ (Lightman 1982; Svensson 1982a, 1984). The dimensionless 
compactness parameter I is defined by I = Lut / RrrieC^ , where L is the luminosity, and R is 
a characteristic length of the object. The above result holds also in a scenario considering 
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injection of high energy particles, which lose their energy via IC scattering of low energy 
photons, and thermalize before annihilating (Svensson 1987; Lightman & Zdziarski 1987). 
The optical depth for scattering by pairs was found in the above analyses to be 0.1 < r± < 
15. However, t± strongly depends on the comoving compactness and sharply increases 
beyond these values when the compactness increases beyond 10'^ (see Lightman & Zdziarski 
1987; Svensson 1987). When the scattering optical depth r± ^ 1 and ^ <^ 1, solving the 
Kompaneets equation gives a cutoff at normalized photon energy hu/rrieC^ — ^/^± (Sunyaev 
& Titarrchuk 1980). 

The results mentioned above are not directly applicable to GRB plasma. The GRB 
plasma is rapidly evolving and steady state can not be assumed. For example: The electron 
distribution function does not reach thermal equilibrium even at low energy (see § 4); Due 
to expansion, the average number of scattering a photon undergoes is ~ t± rather than 
T^; The expansion (and cooling) of plasma electrons has a significant effect on the photon 
spectrum when r± is high. Moreover, significant luminosity at high energies (exceeding the 
pair production threshold) is expected due to synchrotron emission from energetic particles 
in strong magnetic fields, B ~ 10^ — 10^ G, typical for GRB shell-shell collision phase, a 
phenomenon that was not considered in the analyses mentioned above. And last, in the GRB 
case a non-thermal high energy electron population is assumed to be produced, which leads 
at large compactness to the formation of a rapid electro-magnetic cascade, the evolution of 
which was not considered in the past. 

A numeric calculation of GRB spectra that takes into consideration creation and an- 
nihilation of pairs is complicated. The evolution of electromagnetic cascades initiated by 
the annihilation of high energy photons occurs on a very short time scale. On the other 
extreme, evolution of the low-energy, mildly relativistic pairs, which is governed by syn- 
chrotron self-absorption, direct and inverse Compton emission takes much longer time. The 
large difference in characteristic time scales poses a challenge to numeric calculations. For 
this reason, the only numerical approach so far (Pilla & Loeb 1998) was based on using 
a Monte-Garlo method. This model, however, did not consider the parameter space region 
where pairs strongly affect the spectra. Another challenge to numerical modeling is due to 
the fact that at mildly relativistic energies the usual synchrotron emission and IC scattering 
approximations are not valid, and a precise cyclo-synchrotron emission and direct- inverse 
Compton scattering calculations are required. 

In this work, we consider emission within the fireball model framework, resulting from 
internal shocks within an expanding relativistic wind. These shock waves dissipate kinetic en- 
ergy, and accelerate a population of relativistic electrons. We adopt the common assumption 
of a power law energy distribution of the accelerated particles, and calculate the emergent 
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spectra. We present the results of time dependent numerical calculations, considering all 
the relevant physical processes: cyclo-synchrotron emission, synchrotron self absorption, in- 
verse and direct Compton scattering, and pair production and annihilation including the 
evolution of high energy electro-magnetic cascades. A full description of the numerical code 
appears in Pe'er & Waxman (2004). 

We note, that Comptonization by a thermal population of electrons (and possibly 
pairs) was considered as a possible mechanism for GRB production (e.g. Ghisellini & Cclloti 
1999), following the evidence that, at least in some cases, the GRB spectra at low energy is 
steeper at early times than expected for synchrotron emission (Crider et al. 1997; Prcccc et 
al. 2002; Frontcra et a/. 2000). This model is different than the common model considered in 
the previous paragraph in assuming that the kinetic energy dissipated in a collision between 
two "shells" within the expanding wind is continuously distributed among all shell electrons, 
rather than being deposited at any given time into a small fraction of shell electrons that pass 
through the shock wave, and in assuming that energy is equally distributed among electrons, 
rather than following a power-law distribution. There is no known acceleration mechanism 
that leads to the above energy distribution among electrons, and the spectrum predicted 
in this scenario does not account for the claimed steep spectra, which may be naturally 
explained as a contribution to the observed 7-ray radiation of photospheric fireball emission 
(e.g. Meszaros &: Rees 2000). It is nevertheless worthwhile to derive the spectrum that is 
expected from thermal Comptonization under plasma conditions typical to GRB fireballs. 
This will allow to determine whether this process may significantly contribute to GRB 7-ray 
emission. Our numerical code enables an accurate calculation of the pair temperature in this 
scenario, as well as reliable calculation of the emergent spectrum. 

This paper is organized as follows. In §2.1 we derive the plasma parameters at the 
internal shocks stage, and their dependence on uncertain model parameters values. In § 2.2 
it is shown that moderate to large compactness is expected for the parameter range where 
emission peaks at ~ 1 MeV, and approximate analytic results describing the emission in 
this regime are given. Our numerical methods are briefly presented in §3; a detailed de- 
scription of the model is found in Pe'er & Waxman (2004). Our numerical results for the 
scenario of acceleration of particles in shock waves are presented in §4. In §5 we present both 
analytical and numeric calculations of the spectra resulting from thermal Comptonization. 
We summarize and conclude our discussion in §6, with special emphasis on observational 
implications. 
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2. Plasma peirameters and large compactness behavior 

Variability in the Lorentz factor of the relativistic wind emitted by the GRB progenitor 
leads to the formation of shock waves within the expanding wind, at large radii compared 
to the underlying source size. Denoting by F the characteristic wind Lorentz factor and 
assuming variations AF/F ~ 1 on time scale At, shocks develop at radius 2F^cAt = 
5.4 X 10^^ F| 5 Ai_4 cm. For AF/F ~ 1 the shocks are mildly relativistic in the wind frame. 
For our calculations, we consider a collision between two uniform shells of thickness cAt, 
in which two mildly relativistic (F^ — 1 ~ 1 in the wind frame) shocks are formed, one 
propagating forward into the slower shell ahead, and one propagating backward (in the wind 
frame) into the faster shell behind. The comoving shell width, measured in the shell rest 
frame, is AR = TcAt, and the comoving dynamical time, the characteristic time for shock 
crossing and shell expansion measure in the shell rest frame, is tdyn = FAt. 

Under these assumptions, the shocked plasma conditions are determined by 6 model pa- 
rameters. Three are related to the underlying source: The total luminosity L = 10^^ L52 erg s~ 
the Lorentz factor of the shocked plasma F = 10^'^ F2.5, and the variability time At = 
10^'^ At-4 s. Three additional parameters are related to the coUisionless-shock microphysics: 
The fraction of post shock thermal energy carried by electrons ee = 10"*^-^ ee,-o.5 and by 
magnetic field = 10~°'^ e^^-o.s, and the power law index of the accelerated electrons en- 
ergy distribution, d log ^e/d log £e = —P- In the following calculations spherical geometry is 
assumed. However, the results are valid also for a jet like GRB, provided the jet opening 
angle 9 > F~^ (in which case L should be regarded as the isotropic equivalent luminosity). 

In § 2.1 below we derive the characteristic values of plasma parameters obtained under 
the adopted model assumptions, and the characteristic synchrotron and self-absorption fre- 
quencies. Since, as we show in § 2.2, the compactness parameter is related to the optical 
depths for both pair production and scattering by pairs, emission from plasma characterized 
by small compactness is well approximated by the optically thin synchrotron-self-Compton 
emission model. This model, however, ceases to be valid for large value of the compactness. 
In § 2.2 we give an approximate analysis of the emission of radiation from moderate to large 
compactness plasma. 

2.1. Plcisma peirameters 

Denoting by 9p the average proton internal energy (associated with random motion) 
in the shocked plasma, measured in units of the proton's rest mass, the comoving proton 
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density in the shocked plasma is given by 

~ 1 for mildly relativistic shocks and is limited within the fireball model framework to 
Op < few, since the Lorentz factors of the internal shocks can not be larger than a few. This 
is due to the fact that the Lorentz factors of colliding shells can not differ by significantly 
more than an order of magnitude: Shells' Lorentz factors are limited to the range of ~ 100 
to few thousands, where the lower limit of ~ 100 is set by the requirement to avoid too 
large optical depth, and the upper limit of few x 10'^ is due to the fact that shells can not be 
accelerated by the radiation pressure to Lorentz factors ^ 10^ (e.g. §2.3 in Waxman 2003). 

A fraction of the internal energy density Uint — L/ (47rr^cr^) is assumed to be carried 
by the magnetic field, implying 



Equating the particle acceleration time, tacc — ^/icQB) and the synchrotron cool- 
ing time tcooi,sync gives the maximum Lorentz factor of the accelerated electrons, ^ymax — 
{3/2)meC^{q^B)-'^/^ = 6.9 x 10^ L^^^^ e~]L\,^ V^ll At]^^, and the maximum observed energy 
of the synchrotron emitted photons. 



?^7L.T^ = 7 X 10^°(1 + z)-' T,,, eV . (3) 



The minimum Lorentz factor of the power law accelerated electrons, given by 

hAfe)log-(g-) (p = 2) 

is much larger than 7c, the Lorentz factor of electrons that cool on a time scale that is equal 
to the dynamical time scale, which is 7c ~ 1. 

For a typical value of log (£e,moa;/£e,mm) — 7, synchrotron emission from the least ener- 
getic electrons peaks at 

^ J 10^(1 + z)-' <_o.5 6i(%.5 r^i At-_\ el eV (p = 2); 

"-"^ I 5.5 X 10^ (1^)' (1 + L^f <_o.5 ef_,, T^^ Atzl 0^ eV {p > 2). ^> 

The self absorption optical depth t^, — a,^rcAt, calculated using the absorption coefficient 

_ , L3 X lO^^iz-^ Lj^ T^l^ ec,-o.5 €3,-0.5 Ar^cm-^ (p = 2); 

- W.3 X lO-.-V^ L%' r-f/^ e^,_o. ^tzT cm- (p = 3) 
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(Rybicki & Lightman 1979), is smaller than 1 at Speak, 

_ j ^52 <=e.-0.5 ^B,-0.5 ^2.5 ^ (P = 2); , , 

w.. - I ^ ^^1/2^ (p = 3). 

If the fraction of thermal energy carried by the magnetic field is very small, < 
10^^ ^2.5 At_4 e~Q5, the electrons are in the slow cooling regime (i.e., ■jrnin < 7c); the 
power radiated per unit energy below Epeak is proportional to {e/epeakY^^, and the energy 

3/5 

below which the optical depth becomes larger than 1, Sgsa — ^peakTgsa peak 



Ob. 

ssa 



5 X 103(1 + z)-^ Lf 6-io.5 4(^5 ^2^" ^tzl eV {p = 2); 
2.5 X 103(1 + z)-' L'i e-\, ^ Ar^ Qp'''' eV (p = 3). 



When cooling is important (7mm > 7c), as is the case for typical fireball parameters, the 
electron energy distribution, and hence the self- absorption frequency, are modified. As we 
show below (§ 4), for large compactness the energy distribution of electrons and pairs is 
quasi Maxwellian. For a thermal distribution of electrons and positrons, with normalized 
temperature Q = kT/rrieC? and normalized pair density / = n±/np, the self absorption 
frequency vt is approximated by (using the results of Mahadevan, Narayan & Yi (1996) for 
cyclo-synchrotron emission) 

z.T = 5xlO^^ L',^e'^\,V,f"Mz\e^,fl" Hz, (9) 

where 9 = 0.1^_i, and / = lO/i. This result is accurate to better than 10% for parameters in 
the range 0.001 < L52 < 10, 0.01 < < 0.33, 100 < T < 1000,0.1 < ^ < 5, 1 < / < 100. 
The values of 6 and / were found numerically (§ 4) to be within these limits, for a wide 
range of parameters that characterize GRBs. Therefore, 

<t,thermal -huTm + z)-'^ 600(1 + Z)-' L^^O^f^^ AC^ /^/^ eV . (10) 



2.2. Large compactness behavior 

The comoving compactness parameter I' is defined as I' — ARn'^ar, where AR — ctdyn is 
the comoving width and n'^ — e^Lj (47rmec3r^r|) is the comoving number density of photons 
with energy exceeding the electron's rest mass, Sph > nieC^ (in the plasma rest frame). Only 
these photons are of interest, as their number density determines the number density of 
the produced pairs. In deriving the last equation, a mean photon energy (in the comoving 
frame), {Sph) ~ nieC^ is assumed. This assumption is valid as long as the spectral index a 
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{s^dN/de oc £") is not significantly different than 0, and leads to 

^' = IGTrme^c^rSAt = ^50 L52 ee,_o.5 /\tZ\ , . 

— OZU 2 fcpeafe,lMcV;' ^'--2 -L 2.5 ^6,-0.5 ^5,-0.5 >,0.5- 

Here, 6*^ = 10°-^^p,o.5 and At = 10"^At_2 s. Eq. (11) also imphes, using Eq. (5), 

s;lk = 0.3^^ I'T lT ^t-_r e^,, el,, MeV, (12) 

where = 10^/'2. Eq. (12) implies that emission peaking at ~ 1 MeV may be obtained with 
small compactness, /' ~ 1, only for very short variability time, At < lO^"^'^ s, and, using 
Eq. (11), large F, F > 10'^. For the longer variability time commonly assumed in modelling 
GRBs, At ~ 1 to 10 ms (e.g. Piran 2000; Meszaros 2002; Waxman 2003), /' > 1 is obtained 
for the parameter range where synchrotron emission peaks at ~ 1 MeV. The main goal of 
the present analysis is to examine the modification of the spectrum due to the formation of 
pairs in moderate to large compactness GRB plasma. 

The following point should be noted here. At in the range of ~ 1 to 10 ms is commonly 
adopted, since ~ 1 ms variability has been observed in some bursts (Bhat et al. 1992; 
Fishman et al. 1994), and most bursts show variability on ~ 10 ms time scale (Woods & 
Loeb 1995; Walker Schaefer & Fenimore 2000). It should be kept in mind, however, that 
variability on much shorter time scale would not have been possible to resolve experimentally, 
and can not therefore be ruled out. 

Large V implies large optical depth to photon-photon pair production, r^^ I' > 1, and 
also large optical depth to Thomson scattering by pairs, t±. In the absence of pair annihi- 
lation, T± ~ 21' . For Z' ^ 1, T± is expected to be large, implying also that pair annihilation 
is important, since the cross section for pair annihilation is similar to gt {va±{y) ~ cgt 
for sub-relativistic relative velocity v). r-j- may be estimated in this follows. As 

we show in § 4, photons and pairs approach in the case of /' ^ 1 a quasi-thermal distri- 
bution with mildly relativistic effective temperature (or characteristic energy), OrrieC^ with 
^ <^ 1. Under these conditions, the production of pairs via photon annihilation, which for 
I' ^ 1 occurs on a time scale much shorter than the dynamical time and may therefore 
be approximated as the rate of energy production (per unit volume) in > mgC^ photons, 
~ {ceL/ Airr'^r'^c) / {niecHdyn) , is balanced by pair annihilation, the rate of which is given by 
~ n^caT- This implies n± ~ l'^^^ /o'Tctdyn and 

T± « r/2_ (13) 

If synchrotron photons (of energy lower than the pair production threshold ~ mgC^) 
and pairs reach a quasi-thermal distribution via Compton and inverse-Compton scattering 
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interactions, then the pair "effective temperature" may be estimated as follows. The energy 
of low energy photons is increased over a dynamical time by a factor ~ exp(4r±^) (note, that 
due to plasma expansion the number of scatterings is r± rather than r^). The energy Eq of 
the lowest energy photons which reach "thermalization" is therefore given by eQexp{4T±6) ^ 
OnieC^. Assuming that sq > Speak and that the synchrotron spectrum is fiat, e^dN/de oc 
the average energy per photon for synchrotron photons in the energy range Sq < e < me(? 
is ~ £Q\og{meC^ / £q) . Since the number of photons is conserved in Compton scattering 
interactions, and since the number of pairs is much smaller than the number of photons 
(t± l'^^'^), conservation of energy imphes eQ\og{meC^ / £q) ^ OnieC^. Using the relation 
£oexp(4r±^) ^ OnieC^, we therefore find A9t± ^ log[4^r± — log(^)], which implies A9t± ~ 2 
over a wide range of values of t± » 1. The observed effective temperature, V6, is therefore 

(1 + zrVOm^c- ^{1 + zr^^m^c^ = 5^ (^^^^^7^) ^^^^ 

For large optical depth, 1, the plasma expands before photons escape. Assuming that 

the electrons and photons cool "adiabatically" , i.e. that the characteristic energy 9 oc V^^^^ 
where V is the specific volume, the characteristic energy of escaping photons is lower than 
given by Eq. (14) by a factor r^^^^ (since the optical depth drops as V~'^^^). The observed 
characteristic photon energy is therefore 



« (1 + .)-^^ « (1 + .)-^m.c^ = 2^ (^) C'"- MeV. (15) 



In the limit of /' ^ cxo we expect the plasma to reach thermal equilibrium. Assuming 
that the fraction of dissipated kinetic energy carried by electrons is converted to thermal 
radiation, the resulting (blue shifted) radiation temperature is 

rr = 0.1(1 + z)-' iT' L%'' AtzT e:]!o'l MeV. (16) 



3. Numerical calculations 
3.1. Method 

The acceleration of particles in the internal shock waves is accompanied by time de- 
pendent radiative processes, which are coupled to each other. In order to follow the emer- 
gent spectra we developed a time dependent model, solving the kinetic equations describing 
cyclo-synchrotron emission, synchrotron self absorption, Compton scattering (e7 — > 67), pair 
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production {e^e~ — > 77) and annihilation (77 e'^e~). Our model follows the above men- 
tioned phenomena over a wide range of energy scales, including the evolution of the rapid 
electro-magnetic cascade at high energies. 

The calculations are carried out in the comoving frame, assuming homogeneous and 
isotropic distributions of both particles and photons in this frame. Relativistic electrons 
are continuously injected into the plasma, at a constant rate and with a pre-determined 
constant power law index p between 7mm and ^max (see eq. 4), throughout the dynamical 
time, during which the shock waves cross the colliding shells. Above 7max) an exponential 
cutoff is assumed. The magnetic field is assumed to be time independent, given by equation 
2. 

The particle distributions are discrctizcd, spanning a total of 10 decades of energy, 
{iPmin = 10^'^ to jf3„iax = 10^). The photou bins span 14 decades of energy, from Xmin ^ 
Smin/nT'eC^ = 10^*^ to Xmax = ^max/'^e(? = lO''. A fixed time stcp is choscu, typically 
10~^-^ times the dynamical time. Numerical integration, using Cranck-Nickolson second 
order scheme for synchrotron self absorption and first order integration scheme for the other 
processes, is carried out with this fixed time step. At each time step we calculate (i) The 
energy loss time of electrons and positrons at various energies (via cyclo-synchrotron emission 
and inverse Compton scattering taking into account the fact that low energy electrons can 
gain energy via direct Compton scattering); (ii) The annihilation time of electrons and 
positrons; (iii) The annihilation and energy loss time of photons. Electrons, positrons and 
photons for which the energy loss time or annihilation time arc smaller than the fixed time 
step, are assumed to lose all their energy in a single time step, producing secondaries which 
are treated as a source of lower energy particles. The calculation is repeated with a shorter 
time steps, until convergence is reached. 

In the calculations, the exact cross sections for each physical phenomenon, valid at all 
energy range are being used. Calculations of the reaction rate and emergent photon spectrum 
from Compton scattering follow the exact treatment of Jones (1968). Pair production rate, 
and spectrum of the emergant pairs are calculated using the results of Botcher & Schlickeiser 
(1997). Pair annihilation calculations are carried out using the exact cross section first 
derived by Svensson (1982b). The power emitted by an electron with an arbitrary Lorentz 
factor 7 in a magnetic field, is calculated using the cyclo-synchrotron emission pattern (see 
Bekefi 1966; Ginzburg & Syrovatskii 1969; Mahadevan, Narayan & Yi 1996). 

A full description of the physical processes, the kinetic equations solved and the numer- 
ical methods used appear in Pe'er & Waxman (2004). 
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3.2. Adiabatic expansion 

Once the two shocks cross the colhding shells, the dissipation of kinetic energy ceases, 
and the compressed shells expand and cool. The thermal energy carried by protons, elec- 
trons (positrons) and magnetic field decreases and converted back to kinetic energy. If the 
scattering optical depth at the end of the dynamical time is small, photons escape the shells 
and the energy they carry is "lost" from the plasma. If the optical depth is large, then 
photons interact with the expanding electrons and positrons, and this interaction affects 
the emerging spectrum. Since the plasma is collisionless, and particles are coTipled through 
macroscopic electro-magnetic waves, the details of the conversion of thermal to kinetic energy 
are unknown. For relativistic plasma at thermal equilibrium undergoing adiabatic expan- 
sion, the pressure is inversely proportional to V^^^, where V is the volume. We assume that 
the plasma expands in the comoving frame with velocity comparable to the adiabatic sound 
speed, c/\/3, that B oc cx: t~'^ and that electrons and positrons lose energy due to 

expansion at a rate de/e — —dV/3V. 

Since our numerical model is spatially "0-dimensional", we calculate the evolution of the 
photon and particle spectrum in the scenario outlined in the previous paragraph, assuming 
a uniform isotropic particle and photon distribution. This calculation does not take into 
account the energy loss of photons due to the bulk expansion velocity of the electrons, which 
implies that photons are more likely to collide with electrons that move away from (rather 
than towards) them. In order to estimate the effect of this energy loss, we have carried out 
the following calculations. 

We have calculated, using Monte-Carlo technique, the evolution of the momentum of 
a mono-energetic beam of photons emitted at the center of an expanding spherical ball of 
thermal electron plasma, until they escape. The plasma ball was assumed to expand with 
radial velocity v{r) = [r / R{t)]c/\^, where R{t) is the ball radius, and its temperature 
was assumed to decrease as ^ cx R^^ from an initial value of = 0.1. The emergent 
photon spectrum is shown in figure 1 for three different initial photon energies, So/meC^ = 
10~^, 10"^, 10'^, and two initial scattering optical depths, r = 10, 25. These calculations were 
repeated omitting the energy loss of the photons due to the bulk motion of the electrons, 
by assuming v — (while keeping the ball expansion and temperature decrease unchanged) . 
Figure 2 shows the ratio of the average energy of emerging photos with and without inclusion 
of energy loss to bulk motion for several initial optical depths, as a function of the initial 
photon energy. This figure demonstrate that the effect of energy loss to bulk expansion is 
not highly dependent on the initial photon energy, and that it leads to reduction of photon 
energy by a factor ~ 3 for initial optical depths in the range of 10 to 100. 

In the numerical calculations presented in § 4 and § 5 we have corrected for the effect of 
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energy loss due to bulk motion by multiplying the emerging photon energy by the (energy 
dependent) factor inferred from the calculations presented in figure 2. This correction is 
applied for photons of energy exceeding the self absorption energy Egga at the beginning of 
the expansion phase. For photons of energy lower than the self absorption energy igsa at 
the end of the expansion phase, where the optical depth drops to unity, we have applied no 
correction. This is due to the fact that photons at these energies are tightly coupled to the 
electrons, they are continuously emitted and absorbed, and this coupling is the dominant 
factor determining the photons' energy. Note, that the self absorption energy decreases 
during the adiabatic expansion, as the density and the magnetic field decrease. At the 
intermediate energy range, egsa to Egsa, we have applied an interpolated correction factor. 

Figure 3 presents an example of the modification of the spectrum due to bulk expansion. 
Since the fractional energy loss is not strongly energy dependent, the correction we apply 
does not lead to significant distortion of the spectrum. Note, however, that since we do 
not take into account the spread in the energy of emerging photons that initially had the 
same energy, see figure 1, but rather apply a single correction factor appropriate for the 
average energy of emerging photons, the emerging spectrum would in reality be somewhat 
"smeared" compared to our calculation. In particular, the annihilation peak that appears 
at ~ 100 MeV is expected to be "smoothed". 



4. Results 

We have shown in § 2.2 that /' S> 1 is obtained for the parameter range where syn- 
chrotron emission in the fireball model peaks at ~ 1 MeV [see Eq. (12)]. In this section we 
present detailed numerical calculations investigating the emission of radiation at moderate 
to large compactness (for completeness, we present in § 4.1 numerical results also for low 
compactness). We have demonstrated in § 2.2 that for large compactness, the characteristics 
of emitted radiation are determined mainly by with weak dependence on the values of 
other parameters [see, e.g., Eq. (15)]. The results presented in § 4.3 and § 4.2 for particular 
choices of parameter values (e.g. At — 10~^ s and At — 10~^ s) with I' ~ 10^ to 10^, are 
therefore expected to characterize the emission also for other choices of parameters with 
similar values of I'. 
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4.1. Low compactness 

Figure 4 shows spectra obtained for low compactness, /' ^ 3. Synchrotron self ab- 
sorption results in quasi-thermal spectrum at low energies, below Sgsa ~ 100 eV — 1 keV . 
Between Cgsa and ^peafc ~ 10 — 100 keV , the spectral index is softer than expected for syn- 
chrotron emission only: vF^, oc i/" with a ~ 0.3 rather than a ~ 0.5. The reason for this 
deviation is related with the particle population, presented in Figure 5. The self absorp- 
tion phenomenon causes particles to accumulate at low to intermediate energies, forming a 
quasi-Maxwellian distribution with temperature 9 ~ 0.5. Therefore, above this energy, the 
electron spectral index is somewhat softer than the spectral index expected without the self 
absorption phenomenon, [dne/dSe oc instead of drie/dSe oc resulting in a steeper 
slope above Sssa- 

Between Epeak ~ 10 — 100 keV and 10 MeV the spectral slope a ~ —0.3 is harder 
than expected {a = —0.5 for p = 3.0) due to significant inverse Compton emission. The 
combined effects, of relatively soft spectrum at low energies, and cooling of particles by both 
synchrotron emission and Compton scattering lead to the creation of a very soft spectrum 
(a ~ 0.1) near the IC peak at high energies (3 — 30 GeV ). It is therefore concluded, that for 
low compactness, I' ^ 3, and es — Cg, the spectrum expected at all energy bands, between 
~ 100 eV — 30 GeV is flat, with spectral index that varies in the range —0.3 < a < -1-0.3. 

The flattening of the spectrum due to both synchrotron and IC scattering decreases 
the dependence on p. As presented in Figure 6, for p = 2.0 the flux rises slowly at all 
energy bands due to IC scattering, while for p — 2.5 it is nearly constant in the 1 keV - 1 
GeV range. 

The dependence of the spectrum on magnetic field equipartition fraction is shown in 
Figure 7, which demonstrates that comparison of the flux at 1 GeV and 100 keV may allow 
to determine the value of e^. 

4.2. High compactness 

Figure 8 presents results for large compactness. When the compactness is large enough, 
Compton scattering by pairs becomes the dominant emission mechanism. The spectrum 
can not be approximated in this case by the commonly used optically thin synchrotron-self- 
Compton model. As demonstrated in Figure 5, electrons and positrons lose their energy much 
faster than the dynamical time scale, and quasi Maxwellian distribution with an effective 
temperature 9 ~ 0.05 — 0.1 is formed. The energy gain of the low energy electrons by direct 
Compton scattering results in a spectrum steeper than Maxwellian at the low energy end. 
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indicating that a steady state did not develop. As shown in Figure 9, the electron distribution 

approaches a Maxwellian at the end of the adiabatic expansion. The self absorption frequency 
^fsa ~ 3 — 10 keV before the adiabatic expansion (see Figure 3) is well approximated by 
equation (10), valid for thermal distribution of electrons. 

The ratio of pair to proton number density at the end of the dynamical time is / = 
n±/np ~ 10 in the calculations shown in Figure 8, in agreement with the analytical approx- 
imations of § 2.2. This ratio is determined by balance between pair production and pair 
annihilation, and leads to optical depth t± ~ /'^/^ (Figure 10), in accordance with the pre- 
dictions of § 2.2. Pair annihilation phenomenon creates the peak at FmeC^ ~ 10^ r2.5 MeV 
for large compactness. 

Intermediate number of scattering < 20, as is the case for At = 10^*^ s (Figures 
3, 10), leads to moderate value of the Compton y parameter. If the electrons distribution 
is approximated as a Maxwellian with temperature ^ pa 0.1 (see Figure 5), the Compton 
y parameter, y ~ A9t ~ 4 9-iTi is not much higher than 1. In this scenario, the number 
of scatterings is not large enough to create a i/Fi, oc i/^ spectrum, expected when Comp- 
tonization is the dominant emission mechanism (note, that the observed spectral indices at 
1 — 10 keV reported by Crider et al. (1997); Frontera et al. (2000) and Preece et al. 
(2002) are even harder than this value). The resulting spectrum has a slope z/Fj, oc z/° with 
a ^ 0.5 between Egsa ~ 3 — 10 keV and Epeak ~ 10 MeV . For At = 10^^ s, the optical 
depth is r± ~ 60 and the Compton y-paramcter is higher, y ^ 25, resulting in a steeper 
slope. However, even in this scenario the slope is a ~ 0.7 and not the limiting value, a = 1. 

Since the dominant emission mechanism is Compton scattering by the quasi-thermally 
distributed particles, the spectrum is independent on most of the physical parameters related 
with particle acceleration. As is seen in Figure 11, the spectra docs not depend on the power 
law index of the accelerated electrons. The dependence on is weak; smaller Sb gives a 
steeper slope between 1 keV — 10 MeV . 

If /' is large, the flux is rising up to 1 MeV , regardless of the value of p. Therefore, 
observing a decrease in the flux between 30 keV - 1 MeV indicates both p > 2 and /' ^ 3. 

4.3. Intermediate compactness 

For I' in the range of few to few tens, synchrotron emission and Compton scattering 
equally contribute to the observed radiation. Here too, the spectrum cannot be approximated 
by simple analytic formulation. Figure 12 shows several examples of spectra obtained for 
moderate 
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Even though a significant number of pairs are created, / ~ 10 for the three scenarios 
presented in figure 12, the optical depth for scattering is approximately 1, and the self 
absorption frequency is Egsa ~ 1 keV . Although the number of scatterings is not large, 
it is sufficient for increasing the energy at which the spectrum peaks by about a factor of 
3 above the synchrotron model prediction, leading to Speak ~ 500 kcV for Cg and near 
equipartition. Lower values of es lead to higher Speak, while lower values of eg lead to lower 

^peak ■ 

The spectral slope in the 1 keV — 1 McV range is soft a ^ 0.3, instead of the expected 
value a ^ 0.5, similar to the spectra produced for small due to similar reasons. Figure 13 
shows that the exact power law of the accelerated electrons has only a minor influence on 
the observed spectra. 

A significant fiux is expected up to > 1 GeV . Unlike the scenario of very small /', the 
flux decreases above Speak, and no second peak due to IC scattering from high energy electrons 
is formed. This phenomenon is due to pair production, that cuts the flux above the energies 
where such a peak would form, at ~ 1 GeV . Therefore, the main characteristics of spectra 
produced by intermediate values of I' are a moderate increase in the flux in the 1 keV to 
1 MeV energy bands, a moderate decrease in the flux in the 1 MeV to 1 GeV band, and 
a sharp cutoff above this energy. 

5. Quasi thermal Comptonization 

We consider in this section the quasi-thermal Gomptonization scenario proposed by 
Ghisellini & Gelloti (1999). This model is different than the one considered in the previous 
section in that (i) It is assumed that the kinetic energy dissipated in two-shell collision is con- 
tinuously distributed among all shell electrons, rather than being deposited at any given time 
into a small fraction of the shell electrons that pass through the shock wave, and (ii) It is as- 
sumed that energy is equally distributed among electrons, rather than following a power-law 
distribution. In this case, no highly relativistic electrons arc produced, and Comptonization 
is the main mechanism responsible for emission above a few keV: Synchrotron emission is 
self absorbed up to frequency ut, providing the seed photons for Gomptonization creating a 
vF^ oc 1/^ spectrum up to Speak- We derive here the constraints on the peak flux energy Speak 
in this scenario, both analytically and numerically. 

Figure 14 presents numerical results obtained for this scenario. The calculations were 
carried out using a modifled version of the numerical model, in which electrons and positrons 
are forced to follow a Maxwellian distribution, ^(7)^7 oc P^'^e'^'^^^d^. The temperature 9{t) 
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is determined self-consistently by the balance of energy injection and energy loss. The results 
are shown for two representative cases in Figure 14. In both cases, electrons and positrons 
reach a mildly relativistic temperature, 9 ~ 0.3, and the spectrum peaks at Speak ~ 30 MeV . 
The spectral slope below Epeak is a 0.5 {vFi, oc z/") for At = 10^^ s, where t-j. ~ 2, and 
a 1 for At = 10"'^ s, where t± ~ 10. In the former case, the scattering optical depth, 
T± ~ 2, is not large enough to produce the « ~ 1 spectrum expected for r± 3> 1. We 
show below, using analytic analysis, that £peak ^ 30 MeV is a generic result of the thermal- 
Comptonization scenario. 

In an expanding plasma characterized by width A/?, the available time for scattering is 
AR/c, and the time between scattering is {l)/c, where (/) is the mean free path. Therefore, 
the number of scattering, N AR/{1), is proportional to r, and the Compton y parameter 
is given by 

y = 4^T = AeARarUpf = AOTcAtaTUpf. (17) 

The synchrotron spectrum is well approximated by a black body radiation up to frequency 
vt-i given in equation 9. Thus, the observed luminosity is given by 

^ = e^Lsync = ey^mev^e^, (18) 
and the observed emission peaks at 

e^ak = eyhuTT. (19) 

In the following calculation, we assume that e^/ii/^r is not larger than the saturation value 
of £°eafc, 4^rmeC^. If this is not the case, for F = 10^'^, 9 10~^ is needed in order to obtain 
^peak ~ ^ MeV . Since y > few is required, this value of 9 implies a very large optical depth, 
T ~ 10^. Such large value is not obtained for parameter values relevant for GRB fireballs, 
and would lead to strong suppression of emitted radiation. 

Assuming that is given by equation 19, since ^ < 1, in order to get y > few 
photons have to undergo a minimum number of scatterings, i.e. r > few. For a minimum 
value of T = 10°'^To.5, using r = ARaxripf, f is given by 

/ = 4 L^i ri, At.,ro.5. (20) 
Eliminating from equations 18, 19 using equation 9, 6* is given by 

/3_l7f-l/9 7- -2/30 -0.3 _o6. -V3p5/9 

— L.iJ 1^52 -0.5Seafe,lMeV ^ 2.5 > \^^) 

and the Compton y parameter (eq. 17) becomes 

y = 4^r = 5.5f i^S^^^^^l^i^o.sr^rr Ariei,,,Mev"'''- (22) 
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Using the values obtained for ut (eq. 9) and 9 (eq. 21) in equation 19, we obtain 

_ %eofc^ _ ,ot-f-l/18 r-8/15 0.15 -o6. ^/^P^^/^Af 

~huT^~ ^B-0.5Seafe,lMeV ^ 2.5 ^*-4- i^'^) 

The very weak dependence of the right hand side on /, allows eliminating it. Taking /^^^^ ~ 1 
one can take the log of both sides and use equation 22 to express /, 

f-n^^ 7--21/20 27/80 p5 A^9/8 06. 3/8 . s 

Equation 20 then gives 

eia, = 200(1 + z)-' L%''e-^'%i^tzTrlS MeV . (25) 

Since already assumes its maximum value, the only way to meet the observed peak flux 
^peak ~ l MeV in this scenario is by assuming a significantly lower value of total luminosity 
or longer variability time At, both inconsistent with the parameter space region for high 
compactness, as well as with the observations. Note, that the observed peak would be 
obtained at an energy somewhat lower than given by Eq. 25, due to pair production which 
leads to a cutoff at VrUeC^ — 150 r2.5 MeV . 

The above analysis shows that in the "slow heating scenario" the flux cannot peak 
below few hundred MeV . This general result is not changed by the adiabatic expansion, 
during which the optical depth decreases according to T{t) — To[Ai?/ (Ai? + vt)]^, where Tq 
is the optical depth before the expansion and v ~ c/ \/3 is the expansion velocity. The total 
number of scattering increases during the adiabatic expansion by a factor of (1 + \/3) pa 3 
compared to the number of scattering at the dynamical time. This factor enters equations 
20 and 24 with nearly the same power (1 in equation 20 , 9/8 on equation 24), leaving the 
flnal conclusion unchanged. 



6. Summciry &c discussion 

Within the flreball model framework, synchrotron emission peaking at ~ 1 MeV may be 
obtained with small compactness, /' ~ 1, only for very short variability time, At < 10~^'^ s, 
and large F, F > 10^ [see Eqs. (12) and (11), of § 2.2]. For the longer variability time 
commonly assumed in modelling GRBs, At ~ 1 to 10 ms (e.g. Piran 2000; Meszaros 2002; 
Waxman 2003), /' ~ 10^ to V ~ 10^ is obtained for the parameter range where synchrotron 
emission peaks at ~ 1 MeV (For smaller compactness, spectra peak at lower. X-ray, energy). 
This result has two main consequences. First, observed GRB spectra are expected to be 
signiflcantly modifled by the presence of pairs. Second, peak energy » 1 MeV can not be 
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obtained for L ~ lO^^erg/s, since it would imply V » 10^ [see Eqs. (12), (16), (15)], in 
which case most of the radiation will not escape due to large optical depth to Thomson 
scattering by pairs. These conclusions are consistent with the conclusions of Guetta, Spada 
& Waxman (2001), and may provide an explanation for the lack of bursts with peak energy 
> 1 MeV. 

It should be noted at this point that GRB observations do not allow, in most cases, 
to identify variability on ~ 1 ms time scale. Rapid variability, on ~ 1 ms time scale, has 
been observed in some bursts (Bhat et al. 1992; Fishman et al. 1994), and most bursts 
show variability on the shortest resolved time scale, ~ 10 ms (Woods & Locb 1995; Walker 
Schaefer & Fenimore 2000). It should be kept in mind, however, that variability on much 
shorter time scale would not have been possible to resolve experimentally, and can not 
therefore be ruled out. 

We have demonstrated in § 2.2 that for /' in the range of 10^ to 10^ the characteristics of 
emitted radiation are determined mainly by /', with weak dependence on the values of other 
parameters [see, e.g., Eq. (15)]. The peak of the specific luminosity is expected to be close 
to ~ 1 MeV, ~ 1(L52/At _3)V5(r/100)-^ MeV [Eq. (15)], and the spectrum is expected to 
differ significantly from an optically thin synchrotron spectrum. 

These conclusions are consistent with the results of our detailed numerical calculations. 
We have presented numeric results of calculations of prompt GRB spectra within the fireball 
model framework, using a time dependent numerical code which describes cyclo-synchrotron 
emission and absorption, inverse and direct Compton scattering, and pair production and 
annihilation. We have shown that the spectral shape depends mainly on the compactness 
parameter, which is most sensitive to the fireball Lorentz factor F, /' oc F^^. For large 
compactness (small F), /' > 100, the spectra peak at ~ 1 MeV, show steep slopes at lower 
energy, e'^dN/de oc with 0.5 < a < 1, and show a sharp cutoff at ~ 10 MeV (see fig. 8). 
The spectra depend only weakly in this regime on the power-law index p of accelerated 
electrons and on the magnetic field energy fraction (see fig. 11). For small to moderate 
compactness (large F), Z' < 10, spectra extend to > 10 GeV (see figures 4, 12). The spectrum 
at lower energy depends only weakly on p (see fig. 6), but strongly on (see fig. 7). For 
magnetic field close to equipartition, spectra peak at ~ 0.1 MeV for I' ~ 10 and extend to 
higher energy with spectral index a — lO'^^'^'^. 

Since moderate to large compactness is required for a synchrotron peak at ~ 1 MeV, the 
effects of pair-production and direct and inverse- Compton scattering by pairs are predicted 
to be large for observed GRBs. In this case, simple analytic approximations of synchrotron- 
self-Compton emission do not provide an accurate description of the emergent spectra. In 
particular, Compton scattering by the pairs, which are accumulated at intermediate energy 
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9 = e/me(? ~ 0.1 (see figure 9), results in a steep slope, 0.5 < a < 1, in the 1 keV to 
1 MeV band. Still steeper slope, a ~ 3, is obtained at lower energies, below the self- 
absorption frequency determined by the quasi-thermal pair distribution, ef^^ ^ iQO-Sio.s j^^y 
(sec Eq. 10). These steep slopes may account for the steep slopes observed at early times in 
some GRBs (see Frontera et al. 2000; Preece et al. 2002; Ghirlanda, Celloti & Ghisellini 
2003). 

The spectra presented in this paper (§ 4) describe the emission resulting from a single 
collision within the expanding fireball wind, for various choices of model parameters (L, At, F 
etc.). Observed spectra arc expected to be combinations of spectra produced by many single- 
shell collisions, each characterized by different parameters. This is due to the fact that at any 
given time a distant observer is expected to receive radiation from many collisions taking 
place at various locations within the wind. Moreover, observed spectra are inferred from 
measurements where the signal is integrated over time intervals longer than that expected 
for the duration of a single collision. A detailed comparison with observations requires 
a detailed model describing the distribution of single-shell collision parameters within the 
fireball wind (see, e.g., Daigne & Mochkovitch 1998; Panaitescu, Spada,& Meszaros 1999; 
Guetta, Spada & Waxman 2001). The construction and investigation of such models are 
beyond the scope of this manuscript. 

We have shown (§ 5, figure 14) that the quasi-thermal Comptonization scenario (e.g. 
Ghisellini & Celloti 1999), in which kinetic energy dissipated in shocks is continuously 
distributed roughly equally among all electrons, leads to high peak energy, > 30 MeV . This 
scenario may not account therefore for observed GRB spectra. 

Clearly, the most stringent constraints on /', and hence on F, will be provided by mea- 
surements of the spectra at high energy ^0.1 GeV. The fluxes predicted by the model are 
detectable by GLAST^. For large compactness, where emission is strongly suppressed above 
0.1 GeV, the model predicts a steep spectrum at low energy, which is weakly dependent on 
other model parameters (figure 11). For small compactness, where strong emission is ex- 
pected above 10 MeV, the low energy spectrum depends mainly on (fig. 7). Comparison 
of the fiux at ~ 1 GeV and ~ 100 keV (available, e.g. with SWIFT^) may therefore allow 
to determine the fireball magnetic field strength. Stringent constraints on p, the spectral 
index of the energy distribution of accelerated electrons, will be difficult to obtain, due to 
the weak dependence of spectra on this parameter. 



^http://www-glast. stanford.edu 
^ http : / / www . swift . psu. edu 
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Fig. 1. — Monte-Carlo calculations of photon spectra emerging following the injection of 
mono-energetic photons at the center of an adiabatically expanding sphere of electrons with 
initial temperature 9 = kT/nigC^ = 0.1 (see § 3.2). Upper panel: initial scattering optical 
depth T = 10, lower panel: r = 25. Results are shown for several initial photon energy 
(normalized to mgC^): £0 = 10^® (left, thin), Eq = 10~^ (intermediate, bold), Eq = 10^ (right, 
thin). Solid lines show the results of complete calculations, while dashed lines show results of 
calculations where energy loss of photons due to the bulk motion of electrons are neglected. 
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Fig. 2. — The ratio of the average energy of photons emerging from an expanding sphere 
(see Fig. 1 caption) obtained when energy loss to bulk motion is included and excluded. 
Results shown for several initial optical depths: r = 10 (solid), r = 25 (dashed), r = 100 
(dash-dotted). 
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Fig. 3. — The effect of adiabatic expansion on emergent spectra in the case of large com- 
pactness {I' = 250): Dotted: spectrum at the end of the dynamical time (before adiabatic 
expansion); Dashed: spectrum after adiabatic expansion; Solid: after correction for energy 
loss due to plasma bulk motion (see § 3.2). 
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Fig. 4. — Time averaged spectra obtained for low compactness parameter. Results are shown 
for L = 10^2 grg, = = lO^o-^, p = 3 (all cases) and: At = lO^^ g , T = 300 (solid); 
At = 10-3 s , r = 600 (dashed); At = 10"^ s , T = 1000 (dash-dotted). The compactness 
parameter is /' = 2.5, 0.8, 0.6, respectively. Luminosity distance di = 2 x 10^*^ and z = 1 
were assumed. 
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Fig. 5. — Particle distribution at the end of the dynamical time. Thick: At = 10 s,p = 3.0, 
r = 300, /' = 250; Thin: At = 10"^ s,p = 3.0, T = 1000, /' = 0.6. All other parameters 
are the same as in Figure 4. Solid: electron distribution, dash-dotted: positron distribution. 
The dotted lines show Maxwellian distributions at temperatures 6 = kT /m^(? = 0.08 (thick) 
and 6 = 0.5 (thin). 
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Fig. 6. — Dependence of spectra on the power law index p of accelerated electrons, low 
compactness case. Results are shown for At = 10~^ s, F = 300, I' = 2.5, and p = 2.0 (sohd), 
2.5 (dotted), 3.0 (dashed). All other parameters are the same as in Figure 4. Spectra depend 
only weakly on p. 
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Fig. 7. — Dependence on the fraction of thermal energy carried by magnetic field, e^, low 
compactness case. Results are shown for At = 10~^ s, F = 600, I' = 0.8. , and = 0.33 
(solid), tB = 10^^ (dashed), cb = 10~^ (dotted). All other parameters are the same as in 
Figure 4. The ratio of fluxes at 1 GeV and 0.1 MeV is a good indicator for the ratio : Ce- 
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Fig. 8. — Time averaged spectra obtained for high compactness. Results are shown for 
L = 10^2 grg, ee = eB = 10-°•^ p = 3 (all cases) and: At = 10"^ s , T = 300, (sohd); 
At = 10^^ s , r = 300. The compactness parameter is /' = 250, 2500, respectively. The 
scattering optical depth at the end of the dynamical time is 13, 56 respectively. The peaks 
observed at ~ 80 MeV result from pair annihilation. Luminosity distance di = 2 x 10^^ 
and z = 1 were assumed. 
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Fig. 9. — Electrons energy distribution, before and after the adiabatic expansion phase for 
At = 10~^ s,p = 3.0, r = 300, /' = 250 (and all other parameters the same as in Figure 
4). Thick: distribution at the end of the dynamical time; Thin: distribution at the end of 
the adiabatic expansion. The dotted line shows a Maxwellian distribution with temperature 
9 = kT/rrLeC^ ^ 0.08. The particle distribution approaches Maxwellian only at the end of 
the adiabatic expansion phase. 
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Fig. 10. — Energy dependent optical depths for pair production and scattering. Solid: 
At = 10-5 s,p = 3.0, r = 300, /' = 2500; Dashed: At = 10"^ s,p = 3.0, T = 300, /' = 250; 
Dash-dotted: At = 10"^ s,p = 3.0, T = 1000, /' = 0.6. All other parameters are the same 
as in Figure 4. 
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Fig. 11. — Dependence of spectra on the power law index p of accelerated electrons, and on 
the fraction of thermal energy carried by magnetic field, e^, for high compactness. Results 
shown for At = 10"^ s, T = 300, /' = 250 and p = 2.0, = 0.33 (solid); p = 3.0, = 0.33 
(dashed); p = 2.0, = 0.01 (dash-dotted). All other parameters are the same as in Figure 
4. Spectra depend only weakly on p and eg- 
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Fig. 12. — Time averaged flux obtained for intermediate compactness. Solid: At = 10 ^ s 
, r = 300; Dashed: At = IQ-^ s ,T = 400; Dash-dotted: At = 10"^ s , T = 600. All the 
other fireball model parameters are the same as in figure 4. The compactness parameter is 
/' = 25, 6, 8, respectively. 
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Fig. 13. — Dependence of spectra on the power law index p of accelerated electrons, inter- 
mediate compactness. Results shown for At = 10~^ s, F = 300, /' = 25, and p = 2.0 (sohd), 
2.5 (dotted), 3.0 (dashed). All other parameters are the same as in Figure 4. 
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Fig. 14. — Time averaged flux for tfie quasi tliermal Comptonization scenario. Fireball 
model parameters assumed are L = 10^^ erg s~^, F = 300, = = 0.33. At = 10~^ s 
(dashed), At = 10"^ s (solid). 



